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Abstract 

In this paper we present a simple algorithm for representation of 
statistical data of any origin by complex probability amplitudes. Nu- 
merical simulation with Mathematica-6 is performed. The Bloch's 
sphere is used for visualization of results of numerical simulation. On 
the one hand, creation of such a quantum-like (QL) representation 
and its numerical approval is an important step in clarification of ex- 
tremely complicated interrelation between classical and quantum ran- 
domness. On the other hand, it opens new possibilities for application 
the mathematical formalism of QM in other domains of science. 
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1 Introduction 

The problem of interrelation between classical and quantum random- 
ness is extremely complicated, see von Neumann pQ for the pioneer 
study. It has not yet been solved completely, see e.g. Svozil [2], [3] 
for discussions. By the orthodox Copenhagen viewpoint these are two 
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completely different types of randomness. The gap between them is 
huge. However, we can mention a number of attempts to minimize 
this gap, e.g. introduction of Wigner distribution and later develop- 
ment of quantum tomography, see e.g. [5J, [5] and references hereby, 
see also [6] on recent numerical simulation for quantum correlations. 
Mentioned studies were oriented toward classical probabilistic inter- 
pretation of quantum randomness. 

In this paper we would like to formulate and study "inverse Born's 
rule problem" - IBP: to construct a representation of probabilistic 
data by probability amplitudes which match Born's rule. In [7], [8] 
we proposed an algorithm - quantum-like representation algorithm 
(QLRA) - which transforms probabilistic data of any origin in prob- 
ability amplitudes (complex and hyperbolic). In this paper we sim- 
plify QLRA essentially (in particular, by restricting its domain of ap- 
plication). We performed extended numerical simulation based on 
Mathematica-6. The Bloch's sphere is used for visualization of re- 
sults of numerical simulation. On the one hand, creation QLRA and 
its numerical approval is an important step clarification of extremely 
complicated interrelation between classical and quantum randomness. 
On the other hand, it opens new possibilities for application the math- 
ematical formalism of QM in other domains of science, see e.g. [9], 
|12j . By representing probabilistic data by QL-states one might try 
to apply methods of quantum information theory to study such data. 
Of course, only the first step has been done in this direction. 

2 Inversion of Born's rule 

We consider the simplest situation. There are given two dichotomous 
observables a = ct\, ot2 and b = /3\, fa- They can be physical (classical 
or quantum) observables or e.g. two questions which are used for 
tests in psychology, cognitive or social science and so on. We set 
X a = {0:1,02} and A;, = ~ "spectra of observables." 

We suppose that there is given the matrix of transition probabil- 
ities P 6 ' a = (pp a )i where pp a = P(b = 0\a = a) is the probability to 
obtain the result b = (3 under the condition that the result a = a has 
been obtained. 

There are also given probabilities p a a = P{a = a), a £ X a , and 
pp = P(b = /?),/? G Aft. Probabilistic data C = {p%,p b p} is related to 
some experimental context (in physics preparation procedure). It is 
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not assumed that both observables can be measured simultaneously. 
Thus in general two samples (prepared under the same complex of 
conditions) are used to collect the probabilistic data for observations 
of a and b, respectively. 

Our aim is to represent this data by a probability amplitude VO 
(in the simplest case it is complex valued) such that Born's rule holds 
for both observables: 

p b f) = \(Aefy\ 2 ,p a a = \(iP,e a a )\ 2 , (1) 

where {e b ^}/3 & x b and {^a}aex a are orthonormal bases (which are also 
produced by QLRA) for observables b and a, respectively. These ob- 
servables are represented by operators b and a which are diagonal in 
these bases. 



3 QLRA 

In [7], [8] we derived the following formula for interference of proba- 
bilities: 



a 



where the "coefficient of interference" 

2 VUaPaPPa 

This is a trivial mathematical identity. To prove it, one should just 
put A given by (J3]) into ([2]). A similar representation we have for the 
a-probabilities. To simplify considerations, we shall proceed under the 
condition: 

DS: The matrix of transition probabilities P b ' a is doubly stochas- 
tic0 

We also suppose that probabilistic data C = {PaiP'p} consists of 
strictly positive probabilities. We proceed under the basic assumption 
(specifying the type of representation): 

RC: Coefficients of interference Xp,f3 E Xf,, are bounded by one: 

\M < i- 



-'^In a doubly stochastic matrix all entries are nonnegative and all rows and all columns 
sum to 1. 
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Probabilistic data C such that RC holds is called trigonometric, be- 
cause in this case we have the conventional formula of trigonometric 
interference: 

p h p = yi p a ^ a + 2 c ° s ^ J n p >d a - ( 4 ) 

a y a 

where 

A/3 = COS <^g. 

This is simply a new parametrization: a new parameter is used, 
instead of A. Parameters (f>n are said to be b\a-relative phases for the 
data C. We defined these phases purely on the basis of probabilities H 
From the probabilistic viewpoint the formula for interference of 
probabilities is nothing else than perturbation of the classical formula 
of total probability (FTP). We recall this law in the simplest case of 
dichotomous random variables 

p(b = 0) = P(a = ai)P(b = 0\a = a x ) + P(a = a 2 )P(b = (3\a = a 2 ) 

(5) 

Thus the probability P(b = (3) can be reconstructed on the basis of 
conditional probabilities P(b = (3\a = a). FTP plays the fundamental 
role in classical statistics and decision making. However, it is violated 
in experiments with quantum systems. 

We denote the collection of trigonometric probabilistic data by 
the symbol C tr . By using the elementary formula: D = A + B + 
2y/ABcos6 = \\fA + e i(t> VB\ 2 , for real numbers A,B > 0,6 G [0,2vr], 
we can represent the probability p b ^ as the square of the complex am- 
plitude (Born's rule): 

v b P = \m\ 2 - (6) 

Here 

= + VK^T, P G X b . (7) 

The formula ([7|) gives the quantum-like representation algorithm 
- QLRA. For any trigonometric probabilistic data C QLRA produces 
the complex amplitude ip. This algorithm can be used in any domain 
of science to create the QL-representation of probabilistic data. 

2 We have not started with any linear space; in contrast we shall define geometry from 
probability. In the conventional quantum formalism the formula of interference of probabil- 
ities is derived starting directly with the Hilbert space. We recall that in QM interference 
of probabilities is derived via transition from the basis for the a-observable to the basis for 
the 6-observble. From the very beginning observables are given by self-adjoint operators. 
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We denote the space of functions: ip : — > C by the symbol = 
$>(Xb, C). Since X = {/3i,(fo}, the $ is the two dimensional complex 
linear space. By using QLRA we construct the map 

jb\a . C tr _^ 

$(X, C) which maps probabilistic data into complex amplitudes. The 
representation flSJ) of probability is nothing else than the famous Born 
rule. The complex amplitude ip(j3) can be called a wave function of 
the data C or a (pure) state. 

By using the terminology of quantum information theory we can 
say that QLRA represents probabilistic data (of a special sort, namely, 
trigonometric) by qubits. 

We set e|g(-) = £(/? — •)- Dirac delta- functions concentrated in 
points /3 = /?i,/?2- The Born's rule for complex amplitudes © can be 
rewritten in the following form: 

P h p = me h p )\\ 

where the scalar product in the space C) is defined by the stan- 

dard formula: 

<Vi,V2> = 2 MP)MP)- (8) 

The system of functions {e b p}p<zx b is an orthonormal basis in the 
Hilbert space H = ($, (-, •)). 

Let now C R (in general (5 is just a label for a result of observa- 
tion). By using the Hilbert space representation of the Born's rule we 
obtain the Hilbert space representation of the expectation of the ob- 
servable b : Eb = Y,^x b PW)? = Ep e x b PW> 4><^4> = (S^>, 
where the (self-adjoint) operator 6 : H — ► H is determined by its 
eigenvectors: = Peg, (3 G X&. This is the multiplication operator 

in the space of complex functions $(.X&, C) : bip(f3) = Pip(/3). 

To solve IBP completely, we would like to have Born's rule not only 
for the 6-variable, but also for the o-variable: p% = \ (tp, e%)\ 2 , a S X a . 
How can we define the basis {e^} corresponding to the a-observable? 
Such a basis can be found starting with interference of probabilities. 
We have: 

^=VP%Ja 1 + VP%Ja 2 , (9) 

where 

f a _ ( y/Pth^i \ f a _( ^ y/Pfi^ \ nm 
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The condition DS implies that the system of vectors {/„.} is an or- 
thonormal basis iff the probabilistic phases satisfy the constraint: 

4>l3 2 - <A<3i = n mod 27r ' 

Thus, instead of the a-basis (fTUj) which depends on phases, we can 
consider a new a-basis which depends only on the matrix of transition 
probabilities: 

e a ai = ( ) , e" a2 = ( ) (11) 

The a-observable is represented by the operator a which is diagonal 
with eigenvalues ai,«2 in the basis {e^}. 

4 Numerical simulation and visualiza- 
tion on Bloch's sphere 




Figure 1: Transition probability P = 0.1 
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We recall that at the moment we work under condition RC. First 
of all our computer program checks this condition. If RC is violated, 
then the program's output is empty - no point on Bloch's sphere. We 
use spherical coordinates: x = sin 26 cos <j>, y = sin 29 sin <p, z = cos 26. 
A vector on Bloch's sphere is given by 

^ = cos0|O) + sin0e**|l) (12) 

It is convenient to use the QLRA-output in the a-basis. Thus we make 
identification: |0) = e^ x , |1) = e% 2 . We have: p% t = cos#, p% 2 = 

sin#; = cos</>, iyl — Aj^ = svncj). Finally: 
y = ±2 y /p*$- y /l-\* i 

z =P a ai -Pa 2 




Figure 2: Transition probability P = 0.5 



7 



In the progranlj we make the parametrization of probabilities: 
P%i = Q,P% 2 = 1 ~ QiPfa = PiPfr = 1 — P- Bigger values of q give 
more red color, larger values of p give more green. If both q and p 
are rather small, then picture is dark. If both q and p are rather big, 
then picture is yellow. Since we proceed under condition DS, the el- 
ements of the matrix of transition probabilities can be parameterized 

as P/3 iai = P(3 2 a 2 = P, P(3 ia2 = Pfoai = 1 ~ P. 
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